Using high-power lasers for detection of elastic photon-photon scattering 
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The properties of four-wave interaction via the nonlinear quantum vacuum is investigated. The effect of the 
quantum vacuum is to generate photons with new frequencies and wave vectors, due to elastic photon-photon 
scattering. An expression for the number of generated photons is derived and using state-of-the-art laser data 
it is found that the number of photons can reach detectable levels. In particular, the prospect of using the high 
repetition Astra Gemini system at the Rutherford Appleton Laboratory is discussed. The problem of noise 
sources is reviewed, and it is found that the noise level can be reduced well below the signal level. Thus, 
detection of elastic photon-photon scattering may for the first time be achieved. 
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Classically, elastic photon-photon scattering does not take 
place in vacuum. However, according to quantum electrody- 
namics (QED), such a process may occur owing to the interac- 
tion with virtual electron-positron pairs. Several suggestions 
to detect photon-photon scattering have been made, for exam- 
ple using harmonic generation in an inhomogeneous magnetic 
field El, using resonant interaction between eigenmodes in 
microwave cavities 121 |3J, using ultra-intense fields occurring 
in plasma channels j4[, as well as many others, see e.g. Refs. 
JH 01- Related to, but physically different from, elastic 
photon-photon scatterin g is photon splitting JH (see also Jl) 
and Delbriick scatterin g |l(W. the latter being detectable using 
high-Z atomic targets 1 10]. However, no suggestions have yet 
led to detection of elastic scattering among real photons. 

In the present Letter we will investigate the possibility to 
detect photon-photon scattering in vacuum using four-wave 
mixing, where three colliding laser pulses stimulate emission 
in a fourth direction with a new frequency. Four-wave mixing 
has the advantage of not being limited by the low cross sec- 
tion for photon-photon scattering II ill . A similar scheme was 
first studied by R ef. J 1211 . and further theoretical 

El El El 

and experimental II 1611 studies of QED four-wave mixing have 
been performed since then. However, so far the available 
laser technology has not been sufficiently powerful to allow 
for successful detection of scattering events. Moreover, even 
with high laser power, the laser setup and geometry in such 
an experiment is important. In particular, we argue that a 
two-dimensional (2D) setup preferred in earlier literature is 
unlikely to produce a detectable signal. This is in contrast to 
the 3D setup presented in this Letter. We have calculated the 
coupling coefficient for four-wave mixing as a function of in- 
cident angles and laser polarization. For given data of the laser 
beams, together with the coupling coefficients, the number of 
scattered photons can be calculated. The analysis is then used 



to suggest a novel concrete experiment, where the parameters 
are chosen to fit the Astra Gemini (AG) system (operational 
2007) located at the Central Laser Facility (CLF), Rutherford 
Appleton Laboratory (RAL). The estimated number of scat- 
tered photons per shot is 0.07. Furthermore, due to the po- 
larization dependence of the number of scattered photons this 
experimental setup will allow for a unique fingerprint of the 
QED process. We note that a non-perfect vacuum in the in- 
teraction chamber leads to competing scattering events from 
Compton scattering and in principle to collective plasma four- 
wave mixing. Noting that plasma cavitation caused by the 
strong laser pulses will suppress the competing effects rather 
effectively close to the interaction region (see e.g. Ref. fl7ll ). 
we find the noise from competing scattering processes to be 
approximately three orders of magnitude lower than the QED 
signal in our suggested experimental setup. Thus, detection 
of photon-photon scattering will in principle be possible with 
the high-repetition rate AG system. 

An effective field theory containing only the electromag- 
netic fields can be formulated in terms of the Heisenberg- 
Euler Lagrangian El El, which is valid for field strengths 
below the QED critical field 10 16 Vcm -1 and for wavelengths 
larger than the Compton wavelength 10 
Heisenberg-Euler Lagrangian is 
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where S, = he 4 /A5nm 4 c 7 , h = h/2%, h is Planck's constant, e 
is the magnitude of the electron charge, m is the electron mass, 
and c is the speed of light in vacuum. The first part is the 
classical Lagrangian density, while the second part represents 
the lower order non-linear QED-correction. From Q follows 



□E = An [c 2 V ( V • P) - d t (<3,P + cV x M)] 



(2) 
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where □ = d 2 — c 2 V 2 . Here the effective polarization and 
magnetization are given by l2lll 

(3) 
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and 

M={4n)~ 1 % [-2(£ 2 -fi 2 )B + 7(E-B)E] . (4) 

Consider three plane waves, representing the incoming laser 
pulses, with amplitudes allowed to have a weak space time 
dependence, E ; (r,?) = (Ej(r,t)e'^^ +c.c.)/2, where c.c. 
denotes complex conjugate, (pj = ky -r— (Ojt and j = 1,2,3. 
Due to the cubic non-linearity of l|3} and (0), we expect gen- 
eration of a fourth wave with (0)4,154) if the incoming wave 
vectors satisfy the matching conditions 

ki+k2=k3+k4 and 0)1+02 = 0)3 + 0)4, (5) 

for some k|. and a>4 = ck4. 

We assume that the amplitudes change slowly, \dEj/dt | -C 
I Ej j (Oj , j VEj j <C \kj \ \Ej j , so the derivatives in the QED terms 
can be taken to act only on the harmonic parts. Only the res- 
onant terms including the factors E\E\E\ and E^E^Ej, are of 
importance, and hence the wave equation (0 for the generated 
field Eg becomes 

□E g (r,f) =4£o)4 2 G£ 1 £ 2 £3V (k4 ' r ~ fl,4r) (6) 

where G is a geometric factor which only depends on the di- 
rections of the wave vectors and polarization vectors. The 
driving of the initial waves can be neglected due to the weak- 
ness of the generated field, and hence the energy in each of 
these pulses is constant during the interaction. 

The radiation zone solution to the wave equation is 

E g (r,f) = ^Ge i{k ^-^ ,] J / {E X E 2 E^) \ [r e^C*-')*^', 

(7) 

where V' is the interaction region, and the amplitudes are to 
be evaluated at the retarded time tR = t—\r — r l \/c. 

The 2D and 3D geometric factors (G2 C i and G3 f / respec- 
tively) satisfying 0, are given by long and complicated ex- 
pressions 12211 . However, as will be demonstrated below, the 
polarization dependence of Gjd is important for the final de- 
sign of our experiment. 

In practice, the highest power laser systems operates at a 
rather well-defined frequency. However, using frequency- 
doubling crystals we also have access to second harmonics, 
although with some power loss. As will be argued below, this 
makes a 2D geometry less competitive than a 3D setup. As a 
consequence, the 3D configuration of wave vectors which is 
best suited for an experimental setup (allowed by the match- 
ing conditions Q) is given by kj = kx, k 2 = ky, ks = (k/2)z, 
and k4 = kx + ky — (k/2)z. This setup is illustrated in Fig.^ 

Here the wavelength and the direction of the generated 
wave is well separated from that of the incoming laser beams. 
This is highly desirable since optical components will always 
scatter and reflect some light, making it difficult to distinguish 
QED scattered photons of wavelengths close to that of the in- 
coming laser beams. Furthermore, the geometry of the setup 
is not hindered by the fact that the focused laser beams are 
cone-shaped, and the polarizations can be chosen arbitrary 




FIG. 1 : Configuration of the incoming laser beams and the direction 
of the scattered wave for the suggested three-dimensional configura- 
tion of wave vectors, which satisfies the matching conditions {5J. 

since it is not necessary to prevent any counter propagating 
laser beam from entering back into the laser. 

The number of scattered photons can be estimated from 
modeling the laser pulses as rectangular prisms with length 
L and quadratic cross section b , inside which the field ampli- 
tudes are constant during the time of interaction. Although not 
a completely accurate model of the pulse shape, it will give a 
good estimate of the number of generated photons. 

In the 3D case (see Fig.Q the interaction region will take 
the shape of a cube with side b, existing during a time L/c. 
The generated electric field can now be solved for, from which 
the generated intensity and total power can be calculated. Af- 
ter multiplying the generated power by the interaction time 
and dividing by the photon energy T1CO4, the estimated number 
of scattered photons per shot is found to be 

where Gj,d is the geometric factor |G| evaluated for the given 
3D geometry (see Ref. H2I), Pj is the power of the incoming 
pulses, A4 the generated wavelength, and with 77 2 given by 

if = / 112^11 S i n dO d<j) 

J0 JQ Ig, res(n 

where I g is the generated intensity, and 7 gj res its maximum 
value in the resonant direction. The approximate equality of 
T7 2 is accurate within 7% for /7/A4 < 30, i.e. b/X-i < 10 where 
A3 is the fundamental wavelength. 

We note that the scaling of the number of generated pho- 
tons with respect to G^d is decoupled from the pulse model. 
Hence it is possible to test the theoretical predictions by vary- 
ing the polarizations of the incoming pulses, as seen in the 
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upper panel of Fig. [2] where is shown as a function of the 
polarization angle of wave three. For an optimal choice of 
polarization angles, G\ d = 0.77. 

The AG Laser will generate two independently config- 
urable 0.5 PW laser beams of wavelength 800 nm. Each pulse 
will contain a total energy of 15 J, and focused intensities ex- 
ceeding 10 22 WcuT 2 will be reached. The repetition rate is 
expected to be one shot per minute. Using these values, the 
spatial dimensions of the pulse model are taken as b = 1 .6 fim 
and L = 10 jim, which gives tj 2 = 0.025. The 3D config- 
uration described above can be achieved if one of the laser 
beams is frequency doubled and split into two beams. The 
estimated loss of power when frequency doubling is approxi- 
mately 60%, and hence the power of the incoming beams are 
given by Pi = P 2 = 0. 1 PW, P 3 = 0.5 PW. The number of QED 
scattered photons, using optimal polarizations, would then be 
Nm = 0.07 and their wavelength around 267 nm. The lower 
panel of Fig. [3] shows how the scattering number increases 
with increasing laser power, when the focused beam width is 
kept constant. We note that in order to get a statistically suf- 
ficient number of scattered photons, the conditions must be 
reproduced from shot to shot. The pulse energy and length 
fluctuations will be small on a shot to shot basis and vibrations 
are likely to be a minor problem due to the inherent construc- 
tion of the AG system. So, a good shot to shot reproducibility 
is expected from AG, and it can be further improved through 
use of comersially available active control systems. Below we 
will show that already the value = 0.07 for the AG system 
will exceed the estimated noise level. 
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FIG. 2: The upper panel shows the number of scattered photons N^j 
per shot, normalized by the number of photons N max per shot for an 
optimal choice of polarization, as a function of polarization angle <p 
of the wave with direction £3. The lower panel shows N max predicted 
by JsJ when increasing the laser power while keeping the beam width 
constant at b — 1.6 flm, for a system where two source beams are 
used and one of the beams is split into two (solid line) and when three 
source beams are used, hence no beam splitting is required (dashed 
line). 

For comparison, let us consider a 2D configuration where 
two parallel beams collide head on with a third beam, all of 
them with the same wavelength. This configuration gives an 
optimal interaction region and also generates the most scat- 
tered photons, N2d = 2.4 per shot for an optimal choice of 



polarization angles with the AG system. However, the scat- 
tered photons are now emitted along the beam axis and have 
the same wavelength as that of the incoming beams, making 
detection impossible in practice. 

There are other possible experimental configurations that 
we have not investigated in this Letter. It is for instance 
possible to use configurations of two or more lasers of dif- 
ferent wavelengths to generate QED scattered photons with 
wavelengths well separated from the harmonics of the inci- 
dent beams, see Ref. 1 1 6 1 . A 2D configuration consisting of a 
fundamental source beam frequency and its second harmonic 
will not give a scattering direction distinct from the incom- 
ing photon directions. Since present laser systems have suffi- 
cient power only at single frequencies (with second harmon- 
ics available using frequency-doubling crystals) we expect 2D 
configurations to be less competitive. 

A perfect vacuum in the interaction chamber is in practice 
impossible to achieve. Below we give a rough estimate of the 
magnitude of the competing scattering events resulting from 
Compton scattering and collective plasma four-wave mixing, 
in order to determine whether it is possible to distinguish the 
QED scattered photons from noise. 

We assume that we can filter out and detect only the pho- 
tons within the wavelength interval ±50 nm centered at the 
wavelength of the QED scattered photons, 267 nm. We calcu- 
late the scattering cross section for inverse Compton scatter- 
ing into this wavelength interval in order to give a rough esti- 
mate of the noise caused by Compton scattered photons. For 
simplicity, the incident light is assumed to be monoenergetic 
with energy £. We further assume that in the rest frame of the 
electron, the scattering is elastic and isotropic Thomson scat- 
tering, which is a decent approximation as long as ye <§C mc 2 
I23I1 . where 7 is the relativistic factor. 

Using the equations of motion, dp" /dx = —{e/m)F%p b , 
the motion of an electron can be simulated and the average 
cross section for scattering into a given wavelength interval for 
a given field intensity can be estimated, see Fig. [3] Here p a is 
the electron four-momentum, X the electron eigen-time, and 
F a j, is the electromagnetic field tensor. Furthermore, below 
we investigate the effect of the ponderomotive force using a 
Gaussian laserpulse. 

If the electrons were distributed homogeneously within a 
pulse of N p photons, the number of scattered photons during 
a time T would be, N sc = On e cN p X. However, of the scattered 
photons, only a fraction, 8, will be scattered into the detector. 
Electrons accelerated by the laser pulse will most of the time 
have a relativistic velocity in the forward direction of the plane 
spanned by the polarization vector and the wave vector. From 
simulations, and taking into account the relativistic beaming 
effect, we estimate 8 to be approximately 0.01. Of the three 
incoming pulses, only the 800 nm beam will propagate away 
from the detector, see Fig.^ From the above discussion and 
Fig.0 we find that we can neglect the scattering from the 800 
nm beam. The number of photons in one beam is then about 
N p f* 6 x 10 18 (40% of 7.5 J at 400 nm). 

With a maximum intensity of 10 22 WcrrT 2 in the interac- 
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FIG. 3: The average cross section, for inverse Compton scatter- 
ing of incident light with wavelength 800 nm and 400 nm into the 
wavelength interval 216 — 316 nm, as a function of increasing inten- 
sity of the incident light. 



tion region, we find from simulations that the effect of the 
ponderomotive force is crucial. We find that more than 99% 
of the beam photons do not encounter any electrons, due to 
electron blow out, thus we introduce K = 0.01 reducing the 
scattering rate correspondingly. Of the scattered photons, only 
a fraction will have the energies possible to confuse with QED 
scattered photons (within that energy range we introduce the 
subscript A on the cross section a). We find the average cross 
section Oa = 10~ 26 cm 2 from simulations (see Fig. 0. The 
expression for the number of Compton scattered photons per 
beam reaching the detector is then 



N det = K8a A n e cN„T. 



(10) 



We assume the interaction region is imaged onto the detec- 
tor in such a way that the detector only sees a part of the inter- 
action region about 10 jj.m in diameter, which gives an interac- 
tion time of about T~3x 10~ 14 s. We further assume an un- 
perturbed vacuum around 10~ 9 torr, giving n e w 5 x 10 9 cm -3 . 
The number of Compton scattered photons that the detector 
might confuse with QED scattered photons, from both of the 
400 nm beams, is then about A^ e t ~ 5 x 10~ 5 , about three or- 
ders of magnitude lower than the expected number of QED 
scattered photons. Furthermore, because of the high vacuum 
and the electron blow out caused by the ponderomotive force, 
the collective plasma phenomena can be neglected compared 
to the Compton scattering. 

In this Letter, we have investigated the possibility for QED 
elastic photon-photon scattering detection using four-wave 
mixing. The suggested experimental setup consisted of three 
colliding laser pulses in order to stimulate emission of a fourth 
electromagnetic wave with a new propagation direction and a 
new frequency. It has been shown that a 3D configuration, 
in contrast to previously suggested 2D setups, will be able 
to generate a measurable signal for the AG Laser and simi- 
lar high power systems. Moreover, an application of the laser 
polarization makes it possible to obtain a fingerprint charac- 
terizing elastic photon-photon scattering. For given data of 
the laser beams, together with the coupling coefficients, an es- 
timate of the number of scattered photons was obtained. The 
theoretical analysis was then used to suggest a concrete exper- 
iment, where the parameters were chosen to fit the next gener- 



ation AG system at the CLF, RAL. The estimated number of 
scattered photons per shot was 0.07. Due to a non-perfect vac- 
uum in the interaction chamber there is competing scattering 
events such as Compton scattering. However, due to plasma 
cavitation caused by the strong laser pulses the competing ef- 
fects was shown to be suppressed rather effectively. Based on 
ultra-high vacuum technology it was found that the noise from 
competing scattering processes would be approximately three 
orders of magnitude lower than the QED signal in our sug- 
gested experimental setup. Thus, detection of photon-photon 
scattering will in principle be possible with the AG system 
operational in 2007. 
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